Identification of key yeast species and microbe-microbe interactions impacting larval growth of Drosophila in the wild.

Microbiota consisting of various fungi and bacteria have a significant impact on the physiological functions of the host. However, it is unclear which species are essential to this impact and how they affect the host. This study analyzed and isolated microbes from natural food sources of Drosophila larvae, and investigated their functions. Hanseniaspora uvarum is the predominant yeast responsible for larval growth in the earlier stage of fermentation. As fermentation progresses, Acetobacter orientalis emerges as the key bacterium responsible for larval growth, although yeasts and lactic acid bacteria must coexist along with the bacterium to stabilize this host-bacterial association. By providing nutrients to the larvae in an accessible form, the microbiota contributes to the upregulation of various genes that function in larval cell growth and metabolism. Thus, this study elucidates the key microbial species that support animal growth under microbial transition.

Systematic analysis of plasmids of the Serratia marcescens complex using 142 closed genomes.

Plasmids play important roles in bacterial genome diversification. In the Serratia marcescens complex (SMC), a notable contribution of plasmids to genome diversification was also suggested by our recent analysis of >600 draft genomes. As accurate analyses of plasmids in draft genomes are difficult, in this study we analysed 142 closed genomes covering the entire complex, 67 of which were obtained in this study, and identified 132 plasmids (1.9-244.4 kb in length) in 77 strains. While the average numbers of plasmids in clinical and non-clinical strains showed no significant difference, strains belonging to clade 2 (one of the two hospital-adapted lineages) contained more plasmids than the others. Pangenome analysis revealed that of the 28 954 genes identified, 12.8 % were plasmid-specific, and 1.4 % were present in plasmids or chromosomes depending on the strain. In the latter group, while transposon-related genes were most prevalent (31.4 % of the function-predicted genes), genes related to antimicrobial resistance and heavy metal resistance accounted for a notable proportion (22.7 %). Mash distance-based clustering separated the 132 plasmids into 23 clusters and 50 singletons. Most clusters/singletons showed notably different GC contents compared to those of host chromosomes, suggesting their recent or relatively recent appearance in the SMC. Among the 23 clusters, 17 were found in only clinical or only non-clinical strains, suggesting the possible preference of their distribution on the environmental niches of host strains. Regarding the host strain phylogeny, 16 clusters were distributed in two or more clades, suggesting their interclade transmission. Moreover, for many plasmids, highly homologous plasmids were found in other species, indicating the broadness of their potential host ranges, beyond the genus, family, order, class or even phylum level. Importantly, highly homologous plasmids were most frequently found in Klebsiella pneumoniae and other species in the family Enterobacteriaceae, suggesting that this family, particularly K. pneumoniae, is the main source for plasmid exchanges with the SMC. These results highlight the power of closed genome-based analysis in the investigation of plasmids and provide important insights into the nature of plasmids distributed in the SMC.

Global population structure, genomic diversity and carbohydrate fermentation characteristics of clonal complex 119 (CC119), an understudied Shiga toxin-producing E. coli (STEC) lineage including O165:H25 and O172:H25.

Among Shiga toxin (Stx)-producing Escherichia coli (STEC) strains of various serotypes, O157:H7 and five major non-O157 STEC (O26:H11, O111:H8, O103:H2, O121:H19 and O145:H28) can be selectively isolated by using tellurite-containing media. While human infections by O165:H25 STEC strains have been reported worldwide, their detection and isolation are not easy, as they are not resistant to tellurite. Systematic whole-genome sequencing (WGS) analyses have not yet been conducted. Here, we defined O165:H25 strains and their close relatives, including O172:H25 strains, as clonal complex 119 (CC119) and performed a global WGS analysis of the major lineage of CC119, called CC119 sensu stricto (CC119ss), by using 202 CC119ss strains, including 90 strains sequenced in this study. Detailed comparisons of 13 closed genomes, including 7 obtained in this study, and systematic analyses of Stx phage genomes in 50 strains covering the entire CC119ss lineage, were also conducted. These analyses revealed that the Stx2a phage, the locus of enterocyte effacement (LEE) encoding a type III secretion system (T3SS), many prophages encoding T3SS effectors, and the virulence plasmid were acquired by the common ancestor of CC119ss and have been stably maintained in this lineage, while unusual exchanges of Stx1a and Stx2c phages were found at a single integration site. Although the genome sequences of Stx2a phages were highly conserved, CC119ss strains exhibited notable variation in Stx2 production levels. Further analyses revealed the lack of SpLE1-like elements carrying the tellurite resistance genes in CC119ss and defects in rhamnose, sucrose, salicin and dulcitol fermentation. The genetic backgrounds underlying these defects were also clarified.

Dynamic changes in Shiga toxin (Stx) 1 transducing phage throughout the evolution of O26:H11 Stx-producing Escherichia coli.

Shiga toxin (Stx) is the key virulence factor of Stx-producing Escherichia coli (STEC). All known Stxs (Stx1 and Stx2) are encoded by bacteriophages (Stx phages). Although the genetic diversity of Stx phages has frequently been described, systematic analyses of Stx phages in a single STEC lineage are limited. In this study, focusing on the O26:H11 STEC sequence type 21 (ST21) lineage, where the stx1a gene is highly conserved, we analysed the Stx1a phages in 39 strains representative of the entire ST21 lineage and found a high level of variation in Stx1a phage genomes caused by various mechanisms, including replacement by a different Stx1a phage at the same or different locus. The evolutionary timescale of events changing Stx1a phages in ST21 was also determined. Furthermore, by using an Stx1 quantification system developed in this study, we found notable variations in the efficiency of Stx1 production upon prophage induction, which sharply contrasted with the conserved iron regulated Stx1 production. These variations were associated with the Stx1a phage alteration in some cases but not in other cases; thus, Stx1 production in this STEC lineage was determined by differences not only in Stx1 phages but also in host-encoded factors.

Insertion Sequence (IS)-Excision Enhancer (IEE)-Mediated IS Excision from the lacZ Gene Restores the Lactose Utilization Defect of Shiga Toxin-Producing Escherichia coli O121:H19 Strains and Is Responsible for Their Delayed Lactose Utilization Phenotype.

Lactose utilization is one of the general biochemical characteristics of Escherichia coli, and the lac operon is responsible for this phenotype, which can be detected on lactose-containing media, such as MacConkey agar, after 24 h of incubation. However, some Shiga toxin-producing E. coli (STEC) O121:H19 strains exhibit an unusual phenotype called delayed lactose utilization (DLU), in which lactose utilization can be detected after 48 h of cultivation but not after only 24 h of cultivation. Insertion of an insertion sequence (IS), IS600, into the lacZ gene appears to be responsible for the DLU phenotype, and exposure to lactose has been reported to be necessary to observe this phenotype, but the mechanism underlying these phenomena remains to be elucidated. Here, we performed detailed analyses of the lactose utilization abilities of a set of O121:H19 strains and their mutants and found that IS-excision enhancer (IEE)-mediated excision of IS600 reactivates the lacZ gene and that the selective proliferation of IS-cured subclones in lactose-supplemented culture medium is responsible for the expression of the DLU phenotype. In addition, we analyzed the patterns of IS insertion into the lacZ and iee genes in the global O121:H19 population and revealed that while there are O121:H19 strains or lineage/sublineages that contain the IS insertion into iee or intact lacZ and thus do not show the DLU phenotype, most currently circulating O121:H19 strains contain IS600-inserted lacZ and intact iee and thus exhibit this phenotype. IMPORTANCE Insertion sequences (ISs) can modulate gene expression by gene inactivation or activation. While phenotypic changes due to IS insertion/transposition are frequently observed, gene reactivation by precise or simple IS excision rarely occurs. In this study, we show that IS600 is excised from the lacZ gene by IS-excision enhancer (IEE) during the cultivation of Shiga toxin-producing Escherichia coli (STEC) O121:H19 strains that show an unusual phenotype called delayed lactose utilization (DLU). This excision rescued their lactose utilization defect, and the subsequent selective proliferation of IS-cured subclones in lactose-containing medium resulted in the expression of the DLU phenotype. As we also show that most currently circulating O121:H19 strains exhibit this phenotype, this study not only provides information helpful for the isolation and identification of O121:H19 STEC but also offers novel insights into the roles of IS and IEE in the generation of phenotypic variation in bacterial populations.

Global population structure of the Serratia marcescens complex and identification of hospital-adapted lineages in the complex.

Serratia marcescens is an important nosocomial pathogen causing various opportunistic infections, such as urinary tract infections, bacteremia and sometimes even hospital outbreaks. The recent emergence and spread of multidrug-resistant (MDR) strains further pose serious threats to global public health. This bacterium is also ubiquitously found in natural environments, but the genomic differences between clinical and environmental isolates are not clear, including those between S. marcescens and its close relatives. In this study, we performed a large-scale genome analysis of S. marcescens and closely related species (referred to as the 'S. marcescens complex'), including more than 200 clinical and environmental strains newly sequenced here. Our analysis revealed their phylogenetic relationships and complex global population structure, comprising 14 clades, which were defined based on whole-genome average nucleotide identity. Clades 10, 11, 12 and 13 corresponded to S. nematodiphila, S. marcescens sensu stricto, S. ureilytica and S. surfactantfaciens, respectively. Several clades exhibited distinct genome sizes and GC contents and a negative correlation of these genomic parameters was observed in each clade, which was associated with the acquisition of mobile genetic elements (MGEs), but different types of MGEs, plasmids or prophages (and other integrative elements), were found to contribute to the generation of these genomic variations. Importantly, clades 1 and 2 mostly comprised clinical or hospital environment isolates and accumulated a wide range of antimicrobial resistance genes, including various extended-spectrum ß-lactamase and carbapenemase genes, and fluoroquinolone target site mutations, leading to a high proportion of MDR strains. This finding suggests that clades 1 and 2 represent hospital-adapted lineages in the S. marcescens complex although their potential virulence is currently unknown. These data provide an important genomic basis for reconsidering the classification of this group of bacteria and reveal novel insights into their evolution, biology and differential importance in clinical settings.

The global population structure and evolutionary history of the acquisition of major virulence factor-encoding genetic elements in Shiga toxin-producing Escherichia coli O121:H19.

Shiga toxin (Stx)-producing Escherichia coli (STEC) are foodborne pathogens causing serious diseases, such as haemorrhagic colitis and haemolytic uraemic syndrome. Although O157:H7 STEC strains have been the most prevalent, incidences of STEC infections by several other serotypes have recently increased. O121:H19 STEC is one of these major non-O157 STECs, but systematic whole genome sequence (WGS) analyses have not yet been conducted on this STEC. Here, we performed a global WGS analysis of 638 O121:H19 strains, including 143 sequenced in this study, and a detailed comparison of 11 complete genomes, including four obtained in this study. By serotype-wide WGS analysis, we found that O121:H19 strains were divided into four lineages, including major and second major lineages (named L1 and L3, respectively), and that the locus of enterocyte effacement (LEE) encoding a type III secretion system (T3SS) was acquired by the common ancestor of O121:H19. Analyses of 11 complete genomes belonging to L1 or L3 revealed remarkable interlineage differences in the prophage pool and prophage-encoded T3SS effector repertoire, independent acquisition of virulence plasmids by the two lineages, and high conservation in the prophage repertoire, including that for Stx2a phages in lineage L1. Further sequence determination of complete Stx2a phage genomes of 49 strains confirmed that Stx2a phages in lineage L1 are highly conserved short-tailed phages, while those in lineage L3 are long-tailed lambda-like phages with notable genomic diversity, suggesting that an Stx2a phage was acquired by the common ancestor of L1 and has been stably maintained. Consistent with these genomic features of Stx2a phages, most lineage L1 strains produced much higher levels of Stx2a than lineage L3 strains. Altogether, this study provides a global phylogenetic overview of O121:H19 STEC and shows the interlineage genomic differences and the highly conserved genomic features of the major lineage within this serotype of STEC.

Development of a homogeneous time-resolved FRET (HTRF) assay for the quantification of Shiga toxin 2 produced by E. coli.

Shiga toxin-producing Escherichia coli (STEC) is a major intestinal pathogen and causes serious gastrointestinal illness, which includes diarrhea, hemorrhagic colitis, and life-threatening hemolytic uremic syndrome. The major virulence factors of STEC are Shiga toxins (Stx1 and Stx2), which belong to the AB-type toxin family. Among several subtypes of Stx1 and Stx2, the production of Stx2a is thought to be a risk factor for severe STEC infections, but Stx2a production levels vary markedly between STEC strains, even strains with the same serotype. Therefore, quantitative analyses of Stx2 production by STEC strains are important to understand the virulence potential of specific lineages or sublineages. In this study, we developed a novel Stx2 quantification method by utilizing homogeneous time-resolved fluorescence resonance energy transfer (HTRF) technology. To determine suitable "sandwich" assay conditions, we tested 6 combinations of fluorescence-labeled monoclonal antibodies (mAbs) specific to Stx2 and compared the HTRF signal intensities obtained at various incubation times. Through this analysis, we selected the most suitable mAb pair, one recognizing the A subunit and the other recognizing the B subunit, thus together detecting Stx holotoxins. The optimal incubation time was also determined (18 h). Then, we optimized the concentrations of the two mAbs based on the range for linearity. The established HTRF assay detected 0.5 ng/ml of the highly purified recombinant Stx2a and Stx2e proteins and the working range was 1-64 ng/ml for both Stx2a and Stx2e. Through the quantification analysis of Stx proteins in STEC cell lysates, we confirmed that other Stx2 subtypes (Stx2b, Stx2c, Stx2d and Stx2g) can also be quantified at a certain level of accuracy, while this assay system does not detect Stx2f, which is highly divergent in sequence from other Stx2 subtypes, and Stx1. As the HTRF protocol we established is simple, this assay system should prove useful for the quantitative analysis of Stx2 production levels of a large number of STEC strains.

Ddhd1 knockout mouse as a model of locomotive and physiological abnormality in familial spastic paraplegia.

We have previously reported a novel homozygous 4-bp deletion in DDHD1 as the responsible variant for spastic paraplegia type 28 (SPG28; OMIM#609340). The variant causes a frameshift, resulting in a functionally null allele in the patient. DDHD1 encodes phospholipase A1 (PLA1) catalyzing phosphatidylinositol to lysophosphatidylinositol (LPI). To clarify the pathogenic mechanism of SPG28, we established Ddhd1 knockout mice (Ddhd1[-/-]) carrying a 5-bp deletion in Ddhd1, resulting in a premature termination of translation at a position similar to that of the patient. We observed a significant decrease in foot-base angle (FBA) in aged Ddhd1(-/-) (24 months of age) and a significant decrease in LPI 20:4 (sn-2) in Ddhd1(-/-) cerebra (26 months of age). These changes in FBA were not observed in 14 months of age. We also observed significant changes of expression levels of 22 genes in the Ddhd1(-/-) cerebra (26 months of age). Gene Ontology (GO) terms relating to the nervous system and cell-cell communications were significantly enriched. We conclude that the reduced signaling of LPI 20:4 (sn-2) by PLA1 dysfunction is responsible for the locomotive abnormality in SPG28, further suggesting that the reduction of downstream signaling such as GPR55 which is agonized by LPI is involved in the pathogenesis of SPG28.

Genome-wide DNA methylation analysis reveals hypomethylation in the low-CpG promoter regions in lymphoblastoid cell lines.

BACKGROUND: Epidemiological studies of DNA methylation profiles may uncover the molecular mechanisms through which genetic and environmental factors contribute to the risk of multifactorial diseases. There are two types of commonly used DNA bioresources, peripheral blood cells (PBCs) and EBV-transformed lymphoblastoid cell lines (LCLs), which are available for genetic epidemiological studies. Therefore, to extend our knowledge of the difference in DNA methylation status between LCLs and PBCs is important in human population studies that use these DNA sources to elucidate the epigenetic risks for multifactorial diseases. We analyzed the methylation status of the autosomes for 192 and 92 DNA samples that were obtained from PBCs and LCLs, respectively, using a human methylation 450 K array. After excluding SNP-associated methylation sites and low-call sites, 400,240 sites were subjected to analysis using a generalized linear model with cell type, sex, and age as the independent variables. RESULTS: We found that the large proportion of sites showed lower methylation levels in LCLs compared with PBCs, which is consistent with previous reports. We also found that significantly different methylation sites tend to be located on the outside of the CpG island and in a region relatively far from the transcription start site. Additionally, we observed that the methylation change of the sites in the low-CpG promoter region was remarkable. Finally, it was shown that the correlation between the chronological age and ageing-associated methylation sites in ELOVL2 and FHL2 in the LCLs was weaker than that in the PBCs. CONCLUSIONS: The methylation levels of highly methylated sites of the low-CpG-density promoters in PBCs decreased in the LCLs, suggesting that the methylation sites located in low-CpG-density promoters could be sensitive to demethylation in LCLs. Despite being generated from a single cell type, LCLs may not always be a proxy for DNA from PBCs in studies of epigenome-wide analysis attempting to elucidate the role of epigenetic change in disease risks.

CCAT2, a novel noncoding RNA mapping to 8q24, underlies metastatic progression and chromosomal instability in colon cancer.

The functional roles of SNPs within the 8q24 gene desert in the cancer phenotype are not yet well understood. Here, we report that CCAT2, a novel long noncoding RNA transcript (lncRNA) encompassing the rs6983267 SNP, is highly overexpressed in microsatellite-stable colorectal cancer and promotes tumor growth, metastasis, and chromosomal instability. We demonstrate that MYC, miR-17-5p, and miR-20a are up-regulated by CCAT2 through TCF7L2-mediated transcriptional regulation. We further identify the physical interaction between CCAT2 and TCF7L2 resulting in an enhancement of WNT signaling activity. We show that CCAT2 is itself a WNT downstream target, which suggests the existence of a feedback loop. Finally, we demonstrate that the SNP status affects CCAT2 expression and the risk allele G produces more CCAT2 transcript. Our results support a new mechanism of MYC and WNT regulation by the novel lncRNA CCAT2 in colorectal cancer pathogenesis, and provide an alternative explanation of the SNP-conferred cancer risk.